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ABSTRACT 

This paper describes a method of fitting total intensity and polarization profiles in 
VLBI images of astrophysical jets to profiles predicted by a theoretical model. As an 
example, the method is used to fit profiles of the jet in the Active Galactic Nucleus 
Mrk 501 with profiles predicted by a model in which a cylindrical jet of synchrotron 
plasma is threaded by a magnetic field with helical and disordered components. This 
fitting yields model Stokes Q profiles that agree with the observed profiles to within 
the l — 2cr uncertainties; the / model and observed profiles are overall not in such 
good agreement, with the model / profiles being generally more symmetrical than 
the observed profiles. Consistent fitting results arc obtained for profiles derived from 
6 cm VLBI images at two distances from the core, and also for profiles obtained for 
different wavelengths at a single location in the VLBI jet. The most striking success of 
the model is its ability to reproduce the spine-sheath polarization structure observed 
across the jet. Using the derived viewing angle in the jet rest frame, 5' ~ 83°, together 
with a superluminal speed reported in the literature, (3 app = 3.3, yields a solution for 
the viewing angle and velocity of the jet in the observer's frame 6 ~ 15° and (3 ~ 0.96. 
Although these results for Mrk501 must be considered tentative, the combined analysis 
of polarization profiles and apparent component speeds holds promise as a means of 
further elucidating the magnetic field structures and other parameters of parsec-scale 
AGN jets. 
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1 INTRODUCTION 

At radio wavelengths, the jets of active galactic nuclei 
(AGNs) emit synchrotron radiation which is characterised 
by appreciable linear polarization, the plane of polarization 
being perpendicular to the sky projection of the source mag- 
netic field. The polarization structure of jets provides infor- 
mation about the structure of their magnetic fields, which in 
turn influence their evolution and emission properties. Mag- 
netic field structures in AGN jets are of great importance, for 
instance, they affect jet stability. A knowledge of magnetic 
field structure is required in order to translate radio images 
into jet structure and also provides constraints on jet forma- 
tion models. Despite much observational effort the magnetic 
field structures of AGN are not yet well understood. 

Three types of observational results inform our present 
thinking about the magnetic field structures in parsec-scale 
jets. First, there is a tendency for polarization angles to lie 
either parallel or perpendicular to the jet. A tendency for 
BL Lac objects and quasars to display magnetic polarization 
directions respectively perpendicular and parallel to the jets 
has long been known (e.g. Gabuzda et al. 1992; Cawthorne et 



al. 1993), and modern surveys (Lister & Homan 2005) have 
confirmed these general trends. In quasars, there is a broad 
distribution of misalignments but such differences are often 
reduced when Faraday rotation is taken into account (e.g. 
Hutchison et al. 2001). This type of polarization structure 
requires an axisymmetric magnetic field; a helical magnetic 
field is one example, although the production of apparent 
jet magnetic fields that are either parallel or perpendicular 
to the jet direction is not a unique property of helical fields. 

Secondly, transverse Faraday rotation gradients have been 
reported across a number of AGN jets (e.g. Asada et al. 2002, 
Gabuzda et al. 2004). These results, though somewhat con- 
troversial (e.g. Zavala & Taylor 2010), suggest the existence 
of a toroidal magnetic field component. 

Thirdly, a significant number of AGN jets possess obvi- 
ous asymmetries in total intensity and linear polarization or 
other transverse structures that are reminiscent of those re- 
vealed in the helical field simulations of Laing (1981). Note 
that the presence of transverse Faraday rotation gradients 
does not provide any information about whether or not the 
poloidal field component is ordered. Generally, speaking, it 
is only the asymmetry of the transverse intensity and polar- 
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ization profiles that can distinguish observationally between 
a helical field (with an ordered poloidal field component) and 
a toroidal field (with a disordered poloidal field component). 

These observational results lead to the main question ad- 
dressed in this paper: can a helical magnetic field explain 
the observed intensity and polarization profiles of parsec- 
scale AGN jets? 

Of course, asymmetries in the transverse profiles could 
also be attributed to physical asymmetries, such as pres- 
sure gradients or other forms of jet asymmetry. However, a 
helical magnetic field threading the jet of an AGN could po- 
tentially describe the observations without requiring special 
conditions in the jet environment. The main advantage of 
such models is that they can produce the observed trans- 
verse structures while avoiding physical asymmetries which 
might cause the jet to deflect or possibly destabilise. 

In recent years, Laing et al. (e.g. 2006a, 2006b, 2008) have 
also investigated whether ordered magnetic fields are present 
in the jets of FR1 radio galaxies on scales much larger than 
those investigated in this paper. The synchrotron emission 
from various magnetic-field configurations was calculated 
based on three-dimensional models, with the aberration cal- 
culations done via numerical integration for the correct spec- 
tral indices. These theoretical results were compared to the 
observed total intensity and polarized emission from several 
FR1 sources, via fitting to two-dimensional brightness distri- 
butions containing more than 1000 independent pixels. This 
work indicates that the polarization data for larger-scale jets 
are compatible with, but do not require, the presence of an 
ordered toroidal component; at least some kiloparsec-scale 
observations seem to be inconsistent with the presence of 
a strong, ordered, large-scale poloidal component, as is ex- 
pected on theoretical grounds (e.g. Begelman et al. 1984), 
suggesting that large-scale helical fields may be ruled out in 
some cases. However, a transition from a helical field config- 
uration on parsec scales to an ordered toroidal (+ disordered 
poloidal) field configuration on kiloparsec scales is plausible, 
so that a lack of support for the presence of a helical field 
component on kiloparsec scales need not imply the same for 
parsec scales. 

The present paper describes an approach to fitting model 
jet total intensity and polarization profiles to observed pro- 
files taken across VLBI jets. The observed jet profiles for 
the BL Lac object Mrk501 are considered as an example, 
using a model with a cylindrical jet threaded by a heli- 
cal magnetic field, as in Laing (1981). The method pre- 
sented here provides a new means to search for evidence 
for a helical magnetic-field component in parsec-scale jets, 
that does not depend on measurements of Faraday rotation 
gradients, whose reliability can be difficult to determine in 
some cases. The method yields estimates of the viewing an- 
gle and helical-field pitch angle in the jet rest frame, which 
can be used to derive the corresponding quantities in the 
observer's frame when combined with observations of the 
apparent speeds of jet components. 



2 HELICAL MAGNETIC FIELD MODELS 

Laing (1981) investigated three different helical magnetic 
field models, which can be used to predict the total inten- 
sity and polarization profiles across a jet using only two 



parameters, the helical pitch angle 7 and the jet's angle to 
the line of sight 8, both defined in the co-moving frame; 
it is assumed that the velocity of the jet remains constant 
across the jet width. In applying these models to the parsec- 
scale jets of several AGNs, Papageorgiou (2005) found the 
best agreement between observed and model profiles arose 
for the third of these models, in which a helical magnetic 
field of constant pitch angle threads a cylindrical jet. It is 
therefore this model that is the focus of the present work. 

Because the quantity sin 2\ is antisymmetric, where \ is 
the electric vector position angle (EVPA), the contributions 
to the integral of U along the line of sight made by the far 
and near sides of the helical field cancel, so that this integral 
is zero; therefore, the polarization distribution across the 
cylinder corresponds fully to Q (Laing 1981), with 

_ J 90° if Q > 
X ~ \ 0° if Q < 

In other words, the integrated magnetic vector position 
angle is transverse for Q > and longitudinal for Q < 0. 
The derivation of Stokes I and Q for this model can be found 
in Laing (1981). Note that there is a typographical error in 
his final equation for Q(x) in his Appendix A — the sign of 
the second term is incorrect, and the correct expression is 

Q(x) = [a(sin 2 7 — cos 2 7 sin 2 8) — b sin 27 sin 28 
sin 

-csin 2 7(l + cos 2 <S)]. (2.1) 

The analysis of Papageorgiou (2005) showed that this 
helical-field model produced profiles that are considerably 
more strongly polarized than observed. The model was 
therefore modified by introducing a disordered (tangled) 
magnetic field component (see, e.g., Burn 1966). This re- 
quires a third model parameter, the degree of entanglement, 
/, defined as the fraction of the magnetic field energy density 
in tangled form: 

- -J— (2 2) 

(Bl)-l-f (2 ' 2) 

where (Sf.) and (B% } are proportional to the magnetic en- 
ergy densities of the tangled and helical magnetic field com- 
ponents. Increasing the degree of entanglement in the field 
(i.e. increasing /) reduces the degree of asymmetry of the 
total intensity profiles predicted by the model in addition to 
decreasing the degree of polarization. 

Making the same assumptions as Laing (1981) (spectral 
index a = 1, where S oc v~ a , and constant electron density 
throughout the cylinder), profiles for total intensity {I(x)) 
and polarized intensity (Q(x)) across the jet can be derived 
analytically: 

I(x) = \cB% + (1 - f)I L (x) (2.3) 

Q(x) = (l- f)Q L (x) (2.4) 

where Il(x) and Ql(x) are the line-of-sight integrated 
Stokes parameters for Laing's helical-field model, x is dis- 
tance from the jet axis (projected on the sky) and R is the 
radius of the jet cylinder. For Q positive and negative, the 
EVPA is parallel to and normal to the jet direction, respec- 
tively. 
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Figure 1. Transverse structure produced by the model considered, for various viewing angles, 5', and helix pitch angles, 7'. Solid lines 
correspond to total intensity, dark grey regions to longitudinal polarization (EVPAs aligned with the jet) and light grey regions to 
transverse polarization (EVPAs orthogonal to the jet). 



The model described so far describes only the emission 
mechanism in the rest frame of the jet. If the jet is rela- 
tivistic then, although the overall levels of flux density will 
be changed by the Doppler effect and aberration, the trans- 
verse profiles will remain unchanged provided that the jet's 
velocity is constant across its width. 

The Stokes parameters I in the observer's frame (un- 
primed) are related to those in the rest frame of the jet 
(primed) by 

1(8) = l'(8')D a (2.5) 

where 8 is the viewing angle in the observer's frame, 8' is 
the viewing angle in the rest frame of the jet, D = [T(l — 
f3 cos 8)]~ x is the Doppler factor (T is the Lorentz factor and 
/3 the jet velocity divided by the speed of light) and the 
value of a is 2 + a, where a is the spectral index (Scheuer 
and Readhead, 1979; Rybicki and Lightman 1979, Chapter 
4). Strictly speaking, Equation (2.5) refers to the intensity, 
I v , from a continous, optically thin jet, but also applies to 
the flux density per unit length from a cylindrical jet. The 
Stokes parameters Q and U transform in precisely the same 
way. 

If the velocity is constant across the width of the cylin- 
drical jet, the Doppler factor is the same for every point in 
the jet. Therefore, in this case, the Doppler factor will af- 
fect only the overall level of the intensity, not the shapes of 
the I, Q or U profiles. More precisely, the jet profiles I(x), 
Q(x) and U(x) viewed at angle 8 in the observer's frame will 



be identical in shape to the profiles I'(x), Q 1 (x) and U'(x) 
viewed at angle 8' in the jet rest frame. Because the fits were 
obtained by comparing only the shapes of the Q and I pro- 
files, without using the absolute intensities of these profiles, 
the resulting fitted profiles correspond directly to the pitch 
angle and the viewing angle in the rest frame of the jet. Fur- 
ther, the profiles of the polarization angle and the fractional 
polarization will be identical in the two frames. To be ex- 
plicit, a prime is used to denote quantities in the rest frame 
of the jet, and absence of a prime to denote quantities in the 
rest frame of the observer. 

Although the viewing angles in the observer's frame are 
likely to be small, of the order of 1/T, where T is the jet 
Lorentz factor, the rest-frame viewing angles will be much 
larger, and it is these rest-frame values that we have derived 
from our profile fits and quote in the text and tables below. 

This is essentially the approach taken by Claussen-Brown, 
Lyutikov and Kharb (2009), except that they use a force- 
free magnetic field configuration (and a Gaussian jet pro- 
file) defined in the jet's rest frame, whereas a simple helical 
field threading a uniform jet has been used here. This choice 
was made because the objective of this work is precisely to 
look for observational signatures of a helical field compo- 
nent, rather than to look for additional features that will in- 
evitably be present in more sophisticated models. (Note that 
Claussen-Brown et al. (2009) quote observer-frame viewing 
angles, but these can only be given in terms of the (un- 
known) Lorentz factor.) 
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Zakamska, Begelman and Blandford (2008) consider a rel- 
ativistic MHD jet with a toroidal magnetic field. The main 
difference between this model and Laing's is the existence 
of a poloidal field component in Laing's. A purely toroidal 
field configuration can give rise to 'spine-sheath' polariza- 
tion structure and systematic transverse Faraday-rotation 
gradients, but not to asymmetric transverse intensity and 
polarization structure; therefore, it is of interest to try to 
investigate such asymmetric structures using the model pre- 
sented here. 

Broderick and Loeb (2009) and Broderick and McKinney 
(2010) consider theoretical simulations of transverse Fara- 
day rotation measure gradients produced by helical jet mag- 
netic fields. The latter simulations, in particular, directly 
demonstrate the generation of a helical field and associated 
Faraday-rotation gradients due to the combination of the 
rotation of the jet base and the jet outflow. These fully 
relativistic simulations show that the resulting transverse 
Faraday-rotation gradients can sometimes prove to be more 
complex than is predicted by simple helical-field models; 
however, their origin remains the toroidal component of an 
essentially helical jet magnetic field. 

While detailed models for relativistic jets and their mag- 
netic fields have been investigated in these and other studies, 
the least model-dependent way to seek evidence for a heli- 
cal field component is to is to use results for a purely helical 
field; this amounts to focusing on the basic features of trans- 
verse jet profiles rather than aiming to describe them in de- 
tail in the face of many unknown parameters. If high-quality, 
high-resolution data become available for jet profiles, then 
it may become feasible to investigate more detailed models. 



3 MODEL PREDICTIONS 

For convenience, the most notable features predicted by this 
model (Laing 1981), visible in the sample grid of profiles of 
total and polarized intensity for a range of 7' and 5' values 
given in Fig. [I] are summarised here. 

(1) Except for purely toroidal magnetic fields (7' = 90°) 
or viewing angles (in the rest frame of the jet) perpendicular 
to the cylinder axis (<$' = 90°), the distributions of total and 
polarized intensity are usually asymmetric. These asymme- 
tries arise because the change in the sky-projected magnetic- 
field direction occurs most rapidly along the helical field line 
on the side of the jet where the angle between the field and 
the line of sight is smallest. This results in a greater level of 
polarization cancellation along the line of sight on one side 
of the jet than the other. The profile asymmetries are clearly 
seen in cases where neither the helical pitch angle nor the 
viewing angle equals 90° in Fig. [I] 

(2) Displacements between the total intensity profile 
maxima and polarized intensity profile maxima are common. 

(3) The fractional polarization varies considerably across 
the profile. 

(4) The polarized intensity distribution can have one, two 
or three local maxima, and the orientation of the projected 
magnetic field can be either longitudinal, transverse or a 
combination of both within a given profile. 




Figure 2. Effect of finite resolution on the observed magnetic 
field configurations. Convolution was performed using a Gaussian 
beam with a FWHM one quarter the size of the intrinsic jet width. 
Region numbers correspond to the configuration types listed in 
Section 3. 

Further examination of the polarization profiles produced by 
this model shows that 4 different magnetic field distributions 
are possible, as can be seen in Fig. [I] 

(1) Longitudinal all across the jet; e.g., 7' = 10° and 

5' = 30° 

(2) Longitudinal on one side and transverse on the other; 
e.g., 7' = 30° and 5' = 30° 

(3) Longitudinal at the edges and transverse at the cen- 
tre; e.g., i = 30° and 5' = 10° 

(4) Transverse all across; e.g., 7' = 90° and 6' = 90° 

Configuration 2 only occurs when 7' = 5' and configuration 
4 only occurs when 7' = S' = 90° . These special cases would 
not often be expected in nature. However, the effects of finite 
resolution, in which the true jet profile is convolved with a 
Gaussian beam, increase the range of parameter values for 
which configurations 2 and 4 are observed [See Figure 2, 
adapted from Papageorgiou (2005)]. 



4 COMPARISON METHOD 

Papageorgiou (2005) did not carry out any formal model fit- 
ting, and matched the observed and model profiles using a 
number of simple criteria. This approach was time consum- 
ing, and did not necessarily result in the best fit in a 'least- 
squares' sense. This technique has been improved upon here 
by generating a database of theoretical profiles to enable a 
quantitative comparison of the observed and model profiles. 
To generate such a database an estimate of the intrinsic jet 
width must first be made. This was done by generating a 
series of jet profiles of increasing width, which were then 
convolved with the observing beam. The intrinsic jet width 
was then taken to be that for which the convolved I pro- 
file best matched that observed. Both Gaussian and top-hat 
model profiles produced similar results. Since the jet width 
determined in this way is only an estimate, several trial val- 
ues were used, as is described in Section 5. 

The model transverse profiles were generated for a jet 
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Table 1. Best fit Parameters for 6cm Mrk501 Slices (Epoch February 1997) 



Slices 7' 8' f Intrinsic Jet Width 

Slice 1 41° ±3° 81° ±3° 0.70 ±0.05 3.1 mas 

Slice 2 53° ± 1° 80° ± 2° 0.40 ± 0.05 5.7 mas 

Slice 3 47° ± 1° 90° ± 2° 0.00 ±0.10 4.8 mas 




MilliARC SEC 



Figure 3. 6cm Polarization Map of Mrk501 for February 1997 
(Pushkarev et al. 2005). The polarization sticks are proportional 
to the polarized intensity. The three lines across the jet show the 
transverse slices that were analysed. The peak intensity is 0.52 
Jy/beam, and the contour levels are 0.60, 1.25, 2.5, 5.0, 10.0, 
20.0, 40.0, and 80.0% of the peak. 



with angular width determined as described above and con- 
volved with a Gaussian beam corresponding to the observing 
beam, varying the values of 8', 7' and /. The best-fit model 
was taken to be the model giving the smallest residual \ 2 
between the observed and theoretical profiles: 



Q n 



In 



E 



(In ~ 



(4.1) 



where N is the total number of data points, I n and Q n are 
the nth observed total and polarized intensity datapoints re- 
spectively, 7° and Q° are the the nth model total and polar- 
ized intensity model values respectively, and I ma x and Qmax 
are the maximum values of the observed total and polarized 
intensity respectively. The I profiles are downweighted com- 
pared to the Q profiles by the factor Qmax I Imax - This factor 
is close to 0.10 for most of the profiles, and effectively gives 
the Q and I profiles comparable weights in the fitting, while 
also ensuring that the polarization structure of the best fit 
matches the observed polarization structure. In effect, we 
are ignoring the higher signal-to-noise ratios of the I data, 
which is justifiable because the real errors that dominate the 
I image are not due to noise, but instead to mapping errors 
associated with CLEAN, and there is no reason to believe 
these should be much smaller than those for Q. The results 
of the fitting do not depend critically on the specific value 
of this weighting factor. The values of 07 and <tq are given 
by the formulas recently proposed by Hovatta et al. (2012), 
which are based on Monte Carlo simulations, and include 
contributions associated with thermal noise and with uncer- 



tainty introduced by the CLEAN process; the contribution 
of residual D-term uncertainty to <jq is negligible for our 
case, far from the total intensity peak. The thermal-noise 
component was determined in regions far from the region of 
source emission. 

As the number density of electrons and the magnetic field 
strength along the slice are unknown, the database profiles 
were scaled so that the maximum total intensities of the 
observed and model profiles are equal. Matching the to- 
tal flux densities of the observed and calculated I profiles 
would be more accurate; however, doing this would drasti- 
cally increase the computational time required to complete 
the comparisons, as each theoretical profile would have to 
be numerically integrated. In addition, convolution with a 
beam to mimic the effects of finite resolution removes most 
of the asymmetries in the theoretical I profiles. As a re- 
sult, most of the observed I profiles are roughly Gaussian in 
shape. Thus, the scaling factors corresponding to matching 
the observed and theoretical profile maxima or the observed 
and theoretical total fluxes are very similar, justifying our 
approach. 

It is not possible to reliably apply standard statistical 
approaches to estimating the uncertainties in the resulting 
parameters without having well determined uncertainties in 
the fitted quantities — 07 and <jq. Unfortunately, estima- 
tion of the uncertainties associated with values in individual 
pixels of the I and Q images is not straightforward, since 
the imaging process is complex and the values in neighbour- 
ing pixels will be correlated, due to convolution with the 
CLEAN beam. Our 07 and ctq estimates should typically 
correspond to uncertainties in individual pixels that are cor- 
rect to within a factor of a few; nevertheless, these estimates 
may not be good enough to yield fully accurate x 2 values. 
Therefore, although the calculated \ 2 values can certainly 
be used to compare the profiles correspond to yielded by 
different sets of model parameters and identify a set of pa- 
rameters yielding a best fit, the x 2 values cannot be used to 
evaluate the overall goodness of the fits obtained; an alter- 
native method used to obtain estimates of the uncertainties 
in the fitted parameters in Section 5 is described below. 

This method was applied to profiles for slices along which 
the EVPA was either parallel or perpendicular to the jet 
direction. The Stokes parameters Q and U are defined such 
that Q is positive and U — when the EVPA is parallel to 
the jet. 



5 MARKARIAN 501 (MRK 501) 

An ideal VLBI jet for an analysis based on the method de- 
scribed in Section 4 would be one that is straight, well re- 
solved and shows clearly visible transverse I and Q struc- 
ture, with U small compared to Q. Well-resolved VLBI jets 
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Table 2. Best fit Parameters for Mrk 501 Slice 2 at 3 different wavelengths 



Wavelength 


7' 


8' 


/ 


Epoch 


Intrinsic Jet Width 


4 cm 


53° ± 1° 


83° ± 2° 


0.40 ±0.20 


February 1997 


5.5 mas 


6 cm 


53° ± 1° 


80° ± 2° 


0.40 ±0.05 


February 1997 


5.7 mas 


13 cm 


54° ± 1° 


86° ± 2° 


0.35 ±0.05 


May 1998 


21.6 mas 


18 cm 


54° ± 1° 


81° ± 2° 


0.60 ±0.05 


May 1998 


22.9 mas 
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Figure 4. Observed / (solid), Q (dashed) and U (dotted) profiles 
for Slices 1, 2 and 3 (top to bottom) of the 6 cm Mrk501 map 
(see Fig. 3). The scale for the J profile is given along the left- 
hand vertical axis and the scale for the Q and U profiles along 
the right-hand vertical axis. The beam is 1.8 mas along the slice 
and 2.4 mas transverse to the slice. 



are rare, however, and most VLBI jets contain some bends 
(though many of these are most likely small bends amplified 
by projection). Thus, VLBI jets displaying clear transverse 
structure, especially in polarization, were sought, with the 
aim of determining whether profiles across such jets could 
plausibly be represented using the helical field model out- 
lined in Section 2. This requires measuring these profiles 



away from bends, along directions orthogonal to the local 
jet direction. 

The well known active galaxy Mrk 501 (1652+398), which 
has been classified as a BL Lac object with redshift 0.034 
( |de Vaucouleurs et al| [T99T[ > , has been chosen for this initial 
study. The jet of Mrk 501 is almost certainly relativistic, as 
is shown by its one-sidedness; however, the shapes of the ob- 
served jet intensity and polarization profiles should be unaf- 
fected by this relativistic motion. Piner et al. (2010) have 
found apparent component speeds significantly less than 
the speed of light in the Mrk 501 jet. At first sight, this 
would suggest a non-relativistic jet, at variance with the 
one-sidedness of the jet (unless the VLBI jet were intrin- 
sically one-sided). In the context of the standard model for 
VLBI jets, it seems far more likely that these low component 
speeds represent either relativistic motion at a very small an- 
gle to the line of sight or pattern speeds that do not directly 
reflect the speeds of the emitting plasma. In fact, a possi- 
ble detection of superluminal motion with v = (3.3 ± 0.3)c 
based on 43 and 86-GHz VLBI images has been reported 
for this source (Piner et al. 2009). The most detailed multi- 
wavelength studies of Mrk 501 have been carried out by 
Giroletti et al. (2004, 2008). They derived constraints on the 
intrinsic speed of the jet /3 and angle of the jet to the line 
of sight in the observer's frame 8 on various scales, finding 
S < 27° and /3 > 0.88, but with high values /3 > 0.95 allowed 
only for 10° < S < 27°. They also found evidence for limb 
brightening in total intensity on somewhat larger scales than 
those studied here, which they attributed to transverse ve- 
locity structure of the jet. In particular, they proposed that 
the jet has a a fast spine and a slower sheath, as suggested 
by Laing (1996), with the two experiencing different degrees 
of Doppler boosting. The observed limb brightening could 
also have other origins, such as an enhancement in the syn- 
chrotron emission coefficient at the edges of the jet due to 
interaction of outer layers of the jet with the surrounding 
medium, or a helical magnetic field confined to a thin shell 
(e.g. Laing 1981). 

Our selection of Mrk 501 for our analyses is motivated 
by the fact that VLBI images of this AGN show a variety 
of transverse structures that could potentially be associated 
with a helical jet magnetic field, most notably a fairly clear 
'spine-sheath' polarization structure, corresponding to con- 
figuration 3 of Fig. 2 (e.g. Pushkarev et al. 2005). In addition, 
being at a relatively low redshift, the VLBI jet of Mrk 501 
is relatively well resolved, compared to VLBI jets in more 
distant AGNs. 

Fig. 3 shows a 6 cm total intensity map of this source 
with the polarization position-angle sticks superimposed, 
constructed from the same data as those of Pushkarev et 
al. (2005), from February 1997. Near the core, the EVPAs 
are predominantly perpendicular to the jet direction, main- 
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Figure 5. Plots of observed and best-fit model I and Q profiles 
for Slices 1, 2 and 3 (top to bottom) of the 6 cm Mrk501 map 
(see Fig. 3). The observed profiles are solid and the best-fit profiles 
dashed; the gray shaded areas surrounding the observed profiles 
correspond to the range of the lcr uncertainties for Q and the 3(7 
uncertainties for I. The scale for the I profile is given along the 
left-hand vertical axis and the scale for the Q profile along the 
right-hand vertical axis. The beam is 1.8 mas along the slice and 
2.4 mas transverse to the slice. 



taining this orientation as the jet direction begins to change 
about 6 mas from the core. Beyond this region, the EVPAs 
are perpendicular to the jet near the two edges and parallel 
in the centre, forming the spine-sheath polarization struc- 
ture referred to above. The observational results of Giroletti 
et al. (2004, 2008) on comparable scales to those studied 
here are broadly similar to those presented in Fig. 3, though 
at higher resolution, so that the transverse asymmetry is 
more apparent. 

The choice of distances along the jet at which to analyse 
the transverse profiles is constrained by the need to find 
places where the EVPAs are either parallel or perpendicular 
to the jet (i.e., where U is close to close to zero all across the 



Figure 6. Plots of observed and best-fit model I and Q profiles 
for 4 cm (top), 13 cm (middle) and 18 cm (bottom) slices taken 
in the same region as Slice 2 (see Fig. 3). The observed profiles 
are solid and the best-fit profiles dashed; the gray shaded areas 
surrounding the observed profiles correspond to the range of the 
lcr uncertainties for Q and the 3c uncertainties for I. The scale 
for the / profile is given along the left-hand vertical axis and the 
scale for the Q profile along the right-hand vertical axis. The beam 
sizes along and transverse to the slice are 1.4 mas and 1.6 mas 
(4 cm), 3.8 mas and 2.8 mas (13 cm) and 5.1 mas and 3.8 mas 
(18 cm), respectively. 



jet), and the Q profiles have sufficiently high signal-to- noise. 
The profiles to be analyzed were constructed at locations 
where these conditions are satisfied, far from positions where 
the jet bends. A further constraint on choice of position 
arises due to the fact that the jet is essentially unresolved 
in the immediate vicinity of the core. 

The three 6 cm slices chosen for analysis are marked 
on Fig. 3. The profiles were sampled using the AIPS task 
'SLICE', then compared to the database of model profiles 
(each convolved with the Gaussian observing beam) as de- 
scribed in Section 4, and the values of 7', 8' and / resulting 
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in the best fit (minimum residual) identified (given in Table 
1). The intervals between adjacent values of 7', S' and / in 
the database were 1°, 1° and 0.05, respectively. 

The observed / (solid), Q (dashed) and U (dotted) profiles 
for the three 6 cm slices are shown in Fig. 4. The left and 
right sides of the horizontal axis in this figure correspond to 
the North and South sides of the Mrk 501 jet. Recall that 
the Stokes parameters are defined such that Q is positive 
and U is zero for a polarization E field parallel to the jet. U 
must be small (much less than Q) for the model used to be 
valid; the plots show that this condition is satisfied for the 
three slices chosen. 

The observed (solid) and best-fit model (dashed) profiles 
I and Q for the three slices are shown in Fig. 5. Here also, 
the left and right sides of the horizontal axis correspond to 
the North and South sides of the Mrk 501 jet. The I and Q 
curves in Fig. 5 are easily distinguishable, as the Q curves 
are much smaller in amplitude; the intensity scale for I is 
shown to the left and that for Q to the right. The shaded 
bands around the I and Q profiles in Fig. 5 correspond to 3a 
and la uncertainties respectively, estimated in accordance 
with the approach of Hovatta et al. (2012). Unconstrained 
fits for slice 3 yielded best-fit parameters corresponding to 
an opposite sense of the helicity of the magnetic field, com- 
pared to slices 1 and 2; this corresponds to the fact that Q 
changes from negative to positive across slice 2, but becomes 
more negative across slice 3 (Fig. 4). Since such a change in 
helicity between slices 2 and 3 is physically implausible, we 
obtained fits to slice 3 constraining the sense of the helicity 
of the magnetic field to be the same as for slices 1 and 2; 
this is the fit shown in Fig. 5 (and given in Table 1). The 
best-fit model Q profiles for slices 1 and 2 lie within la of the 
observed Q profile, while the best-fit Q profile for slice 3 lies 
within 2a of the observed profile. The model and observed I 
profiles for slices 1 and 2 differ by more than 3a over sustan- 
tial fractions of these profiles, with the model profile being 
more symmetrical than the observed profile; the fitted and 
observed I profiles for slice 3 agree to within about 2a. 

Since the intrinsic (or unconvolved) angular jet width was 
estimated using the procedure outlined in Section 4, the 
fitting procedure was repeated for a series of intrinsic jet 
widths centered on the values used above and listed in Ta- 
ble 1. Changing the intrinsic jet width had only a very minor 
effect on the best fit 7' and 5' values while having a more 
pronounced effect on the degree of entanglement. However, 
in general, as the intrinsic jet width varied from its esti- 
mated value, the x 2 value of the best fit profiles increased. 
While the intrinsic jet width that minimised the % 2 value for 
a given slice is not exactly that found by fitting the / profile, 
the two never differed by more then 10%. Values differing 
by more than 15% yielded very large values of \ ■ 

As was noted above, the uncertainties in the fitted pa- 
rameters cannot be estimated reliably using standard sta- 
tistical techniques, since the values of ai and oq used are 
only estimated uncertainties. Instead, estimates of the pa- 
rameter uncertainties were derived by carrying out Monte 
Carlo simulations, as follows. Infinite resolution model maps 
corresponding to model parameters corresponding to those 
inferred for the Mrk501 jet were generated, and the Fourier 
transforms of these maps sampled at the set of baselines 
corresponding to the observations. This effectively yielded 
model visibilities at the set of UV baselines that were actu- 



ally observed. Noise consistent with the noise levels observed 
in the images was then added. The resulting 'noisy' model 
data were imaged in AIPS in the usual way, profiles were 
taken across the resulting maps and these profiles were sub- 
ject to the same fitting process as the real observed profiles. 
The resulting best-fit parameters were then compared to the 
parameters used to generate the model maps to estimate the 
uncertainty in the fitting process. This was carried out for 
200 Monte Carlo realizations of the model data, and the 
distribution of derived uncertainties used to estimate the 
corresponding la uncertainties, listed in Tables 1 and 2. 

The line-of-sight angles, 8' , for the first two slices in Ta- 
ble 1 agree to within a degree, and differ by less than la, 
as expected if the intrinsic bends in this region of the jet 
are small. The value of 5' for slice 3 differs by about 10° 
(5' = 90°), which corresponds to about 3.5a; however, as 
was pointed out above, it was necessary to constrain the fit 
for slice 3 to have the same sense of helicity as slices 1 and 
2. Therefore, we do not feel confident that the change in 
viewing angle implied by the nominal (constrained) best fit 
for slice 3 is significant. 

The fitting results indicate a somewhat higher value of 7' 
for slice 2 than for slice 1, with this difference being about 
4a, suggesting that the variation in 7' between the slices 
may be real. In this case, this suggests that the appearance 
of a spine-sheath polarization structure at the position of 
slice 2 could be due to a (relatively small) change in the 
helical pitch angle. The value of 7' for slice 3 lies between 
the values for the other two slices; formally, 7' for slice 3 
differs from the value for slice 2 by about 3.5cr, and from 
the value for slice 1 by slightly less than 2a. 

The values of / appear to decline with distance from the 
core, from 0.7 at slice 1 to 0.4 at slice 2. This trend may 
continue for slice 3 (/ = 0), although this is unclear due to 
the uncertainty associated with slice 3 discussed above. The 
difference between slices 1 and 2 appears to be significant: 
the two / values differ by approximately 4a; physically, this 
represents a decrease in the disordered component of the 
magnetic field with distance from the core in the region of 
the jet sampled. 

The fitting procedure was repeated for slice 2 using images 
at wavelengths of 4 cm for February 1997 (data of Pushkarev 
et al. 2005) and 13 cm and 18 cm for May 1998 (data of Crokc 
ct al. 2010). The results are given in Table 2 and shown 
in Fig. 6. It is clear that the observations at these other 
wavelengths yield consistent best-fit parameter values; the 
values of 7', S' and / essentially all agree to within la, with 
the largest differences not exceeding about 2a. 

Fits were also obtained for Slices 1, 2 and 3 using Model B 
of Laing (1981), which consists of a cylinder of fixed radius 
containing a randomly orientated magnetic field with no ra- 
dial component. The only parameter in this model is the 
viewing angle. This model can produce a region of trans- 
verse field at the jet axis surrounded by regions of longi- 
tudinal field, but the profiles produced by this model are 
completely symmetric. Because Model B also produced pro- 
files that are considerably more strongly polarized than is 
observed, we included the degree of entanglement / in these 
models as well. Figure 7 shows the best-fit profiles for the 
three slices for the 6cm image of Mrk501. The model fits for 
slices 1 and 2 are clearly considerably poorer than those in 
Fig. 5. In contrast to Fig. 5, where the model and observed 
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Table 3. Best fit Parameters for 13 cm Mrk501 Slice for three 
values of a 
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Figure 7. Plots of observed and best-fit model I and Q profiles 
for 6 cm slices 1, 2 and 3 obtained for Model B of Laing (i981). 
The observed profiles arc solid and the best-fit profiles dashed; the 
gray shaded areas surrounding the observed profiles correspond to 
the range of the la uncertainties for Q and the 3<r uncertainties for 
I. The scale for the / profile is given along the left-hand vertical 
axis and the scale for the Q profile along the right-hand vertical 
axis. The beam is 1.8 mas along the slice and 2.4 mas transverse 
to the slice. 
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6 SPECTRAL INDEX 

A spectral index (a) of unity was assumed throughout the 
jet in order to make the equations that govern this model 
analytically solvable. However, observations of the spectral 
indices in the jets of AGN have shown that this value is typ- 
ically ~ 0.5. This gives rise to concerns about the trustwor- 
thiness of the fits obtained assuming a = 1.0. To investigate 
this, model profiles were obtained numerically for two other 
spectral indices. 

Profiles for both I and Q were determined numerically by 
integrating equations Al and A2 from Laing (1981) using 

This 



1964) 



the Gauss-Kronrod quadrature (Kronrod, A.S. 
method was required as the functions governing both I and 
Q have singularities at one of their limits. 

Figure[8]shows the effect of variations in the spectral index 
on the I and Q profiles. The left and right hand plots show 
profiles for jets with a = 1 and a — 0.5 respectively, and 
the same values of 7, S and /. The differences between these 
profiles are very small. The positions where the polarization 
angle changes by 90° are slightly different. In addition, the 
amount of longitudinal polarization has increased slightly, 
while the amount of transverse polarization has decreased. 
This second difference results in the configuration map for 
a — 0.5 being slightly different to that for a = 1 (compare 
Fig. 2 and 9). 

The effect of changing spectral index on the best fit model 
values of 7', 6' and / was examined using slices through the 
13 cm images at the location of Slice 2 in Fig. 3. Model 
databases were then generated numerically using values of 
a — 0.507 (the average value of a across the slice) and 
a = 0.248 (the minimum value of alpha across the slice, 
representing the maximum deviation from unity) and the 
FWHM of the beam used in creating the 13 cm map. These 
databases were then used to obtain new best-fit parameters 
for the 13 cm slice. The results are shown in Table 3. The dif- 
ferences in the best fit parameters are extremely small with 
the 7',$' and / values for the two numerically calculated 
best-fit profiles (a = .508 and a = .248) and the analyti- 
cally calculated best-fit profile (for a — 1.0) all coinciding 
within the la errors. 



Q profiles agree to within la everywhere, the differences be- 
tween the model and observed Q profiles in Fig. 7 exceed 
la essentially everywhere, with the differences exceeding 3a 
along most of the profiles. The fit for Model B for slice 3 
is very similar to the corresponding fit for the helical-field 
model in Fig. 5. Thus, while Model B of Laing (1981) can 
provide a comparably good fit for slice 3, it does not provide 
a viable alternative to the polarization structure across the 
other two slices. 



7 COMPARISON WITH RESULTS FROM 
FARADAY ROTATION PROFILES 

As described earlier, the fractional polarization from the he- 
lical field model is generally lower on the side of the jet 
where the field lines are closest to the line of sight. If the 
magnetic field threading the Faraday rotating medium is co- 
moving with the jet and is threaded by the same magnetic 
field, the magnitude of the Faraday rotation measure will 
be highest on the side of the jet where the field lines are 
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Figure 8. Theoretical profiles generated numerically for -y' = 60°, 
a = 1.0 and a = 0.5, respectively. The magnetic-field orientation is 

closest to the line of sight, i.e., on the side with the lower 
fractional polarization. It is therefore of interest to check 
whether the transverse polarization and Faraday-rotation 
profiles observed for the Mrk 501 jet are consistent with 
this picture. In particular, Gabuzda et al. (2004) reported 
the detection of a transverse Faraday rotation measure (RM) 
gradient across the Mrk 501 jet, based on 2+4+6 cm VLBA 
polarization observations, and interpreted this as evidence 
for a helical magnetic field associated with this jet. Croke et 
al. (2010) subsequently reported a transverse RM gradient 
in the same sense based on 4+6+13+18 cm VLBA polariza- 
tion observations. 

Taylor & Zavala (2010) have recently questioned the va- 
lidity of reported transverse RM gradients, and proposed a 
number of criteria for the reliable detection of such gradi- 
ents, namely that there be (i) at least three "resolution ele- 
ments" (taken to mean three beamwidths) across the jet, (ii) 
a change in the RM of at least three times the typical error 
a, (iii) an optically thin synchrotron spectrum at the loca- 
tion of the gradient, and (iv) a monotonic, smooth (within 
the errors) change in the RM from side to side. The more 
recent Monte Carlo simulations of Hovatta et al. (2012) sug- 
gest that a more reasonable limit for the distance spanned by 
the RM values comprising the gradient is two beamwidths 
if the RM values differ by at least 3<j, and as little as 1.5 
beamwidths in the case of higher signal-to-noise ratios. 

As part of their study, Taylor &: Zavala (2010) con- 
sider several RM maps based on multi-wavelength 2-4 cm 
VLBA data, searching for RM gradients specifically in re- 
gions where the RM distributions span more than three 
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8' = 40° and / = 0. The left and right profiles were determined for 
perpendicular to the polarization orientation. 

beamwidths. On this basis, they report an absence of sig- 
nificant transverse RM gradients in the Mrk 501 jet, based 
on profiles taken roughly 10 mas from the core, and therefore 
call into question the earlier results of Gabuzda et al. (2004) 
and Croke et al. (2010). We do not wish to analyze this dis- 
crepancy between the results of Gabuzda et al. (2004) and 
Croke et al. (2010), on the one hand, and Taylor & Zavala 
(2010), on the other hand, in full detail here. However, for 
transparency, we will make some remarks concerning this, 
as a justification for our further consideration of the results 
of Gabuzda et al. (2004) below. 

The transverse RM gradient reported by Gabuzda et al. 
(2004) essentially satisfies the last three criteria of Taylor 
& Zavala (2010), but not the first; however, because the 
gradient spans approximately 1.7 beamwidths, it marginally 
satisfies the modified criterion of Hovatta et al. (2012). The 
gradient reported by Croke et al. (2010) spans roughly four 
beamwidths across the jet, and satisfies all four criteria of 
Taylor & Zavala (2010). Note that the profile considered 
by Taylor & Zavala (2010) is further from the core than 
the region considered by Gabuzda et al. (2004), where the 
magnitude of the RM values would be expected to be lower; 
the region considered by Croke et al. (2010) overlaps with 
the region of the profile of Taylor & Zavala (2010), but the 
data of the former study were considerably more sensitive 
to small amounts of Faraday rotation, since they included 
observations at the longer wavelengths of 13 and 18 cm. 

Thus, while the resolution of the images of Gabuzda et 
al. (2004), analyzed in the same general region as our own 
profiles, may not represent conclusive proof of a uniform 
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Figure 9. Magnetic-field configuration map analogous to Fig. 2 
for a = 0.5 after convolution with a Gaussian beam with a 
FWHM one-quarter the size of the profile to mimic the effects 
of finite resolution. The region numbers correspond to the config- 
uration types listed in Section 3. 
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Figure 10. 6 cm total intensity contours given in Fig. 3 super- 
posed with a colour image of the distribution of the degree of 
polarization. The increase in the degree of polarization toward 
the edges of the jet, also reported by Pushkarev et al. (2005), is 
clearly visible. The locations of the highest degrees of polariza- 
tion do not coincide with those expected for curvature-induced 
polarization. 



transverse RM gradient, their results can reasonably be con- 
sidered at least tentative evidence for a Faraday-rotation 
gradient across the Mrk 501 jet, particularly in light of the 
subsequent, better-resolved results of Croke et al. (2010). 
Therefore, we feel justified in using these results as a work- 
ing hypothesis. We emphasize that we are interested here 
only in the very crude question of whether there seems to 
be a higher typical RM magnitude on one or the other side 
of the VLBI jet, and whether this is consistent with the side 
predicted by our simple model. 

In the regions of the jet considered here, comparing the 
two sides of the jet, the Northern side (corresponding to the 
left-hand sides of the plots in Figs. 5-6) generally displays 
higher fractional polarization, i.e., more negative values of 



Q/I. It follows that the Faraday rotation measures on the 
Northern side of the jet, RMjv, are expected to be smaller 
in magnitude than those on the Southern side, RMs. This is 
consistent with the observational results of Gabuzda et al. 
(2004), based on the same 2-6 cm VLB A data considered 
here: RMjv = -55 ± 20 rad/m 2 , RM S = 130 ± 20 rad/m 2 . 
The effect of the measured Galactic rotation measure in the 
direction toward Mrk 501 [+42 rad/m 2 , Rusk (1988)] was 
subtracted from the observed polarization angles before the 
Faraday rotation map was constructed, so that the residual 
observed rotation measure should correspond to Faraday ro- 
tation occurring in the vicinity of the AGN. If the Galactic 
rotation measure is instead taken to be the value typical of 
the region within a few degrees of Mrk 501 in the catalog of 
Taylor, Stil and Sunstrum (2009), ~ +20 rad/m 2 , these two 
values become RMjv — —35 rad/m 2 , RMs ~ 150 rad/m 2 . 

We can crudely estimate the expected quantitative differ- 
ence between the magnitudes of the Faraday rotation mea- 
sures on either side of the jet if the Faraday rotating material 
is concentrated in a relatively thin shell in outer layers of the 
jet using Eq. (2) of Laing (1981): for 7' = 40° and 8' = 80, 
this indicates a difference of about 35%, somewhat smaller 
than the observed value of about 60 — 75%. 

Considering this one case, and given the limited trans- 
verse resolution of the rotation-measure image of Gabuzda 
et al. (2004), it could be a coincidence that the model con- 
sidered correctly predicts the side of the jet that should 
have the higher Faraday rotation. However, our analysis 
here illustrates the sort of comparison that could, in princi- 
ple, usefully be carried out once profile-fitting results for a 
greater number of AGNs are available. Of course, the sign 
of the Faraday rotation measures cannot be predicted by 
our model, as the polarization profiles do not depend on the 
polarity of the magnetic field. 



8 DERIVATION OF 8 AND p = V/C IN THE 
OBSERVER'S FRAME 

The line of sight angle 8' in the rest frame of the jet is related 
to the corresponding angle in the observer's frame, 8, by 



sin 8' 



sin 8 



T[l - p cos 8] 



(8.1) 



where F is the Lorentz factor for the motion and P = v/c 
(Rindler 1990). 

In addition, the apparent speed of a component moving 
along the jet can be written 

P sin 8 



1 — p cos 8 
TP sin 8' 

{p- 2 - ir 1/2 sin<5' 



(8.2) 

(8.3) 
(8.4) 
(8.5) 



It follows that, if p apP is known from component motion 
and 8' is known, for example, from profile fitting such as that 
carried out in this paper, Equation (8.4) can be rearranged 
to give 



1 + 



• 2 s n -V2 
sin 8 
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The profile fitting has indicated 8' ~ 83° for the Mrk 501 
jet. Taking this together with the superluminal motion re- 
ported by Piner et al. (2009), /3 app = 3.3, gives /3 ~ 0.96 
and 8 ~ 15°. Using the range of 8' values given by the var- 
ious profile fits, 80° — 86°, does not change the resulting 
values of 8 by more than 1°. These results are consistent 
with the conclusion of Giroletti et al. (2004), based on com- 
pletely different information, that 8 < 27° and /3 > 0.88, 
with > 0.95 allowed only for 10° < 8 < 27°. Although 
these particular results are somewhat uncertain, since the 
relevant superluminal motion corresponds to a weak feature 
whose position is difficult to determine precisely, this illus- 
trates the potential of this approach. Note as well that using 
the more typical subluminal speeds obtained for the jet of 
Mrk 501, e.g. /3 app ~ 0.47 (Piner et al. 2010), implies angles 
to the line of sight 8 ~ 60° and /3 ~ /3 app , which are not rea- 
sonable, since they cannot provide the observed one-sided 
VLBI structure. This implies that the very low component 
speeds typically observed in Mrk 501 are incompatible with 
its observed polarization structure, and so must represent 
pattern speeds rather than physical speeds (in other words, 
they do not represent highly relativistic motion viewed at a 
very small angle). 

The approach taken in this section is similar to that used 
by Canvin et al. (2005) on kiloparsec scales, in that both 
analyses fit a jet model to the observed profile in order to 
determine the angle between the jet axis and the line of sight 
in the rest frame of the flow. One advantage of applying this 
technique on parsec scales is that it is possible to obtain 
constraints on the underlying flow speed fairly directly using 
observations of superluminal components, whereas, in the 
kiloparsec case, more indirect estimate must be used, such 
as the jet-counterjet brightness ratio. 



9 DISCUSSION 

In general, the procedure described above has yielded good 
fits to the polarization data for Mrk 501. The deviations 
between the model and observed Q values for slices 1 and 
2 are less than the estimated la uncertainties for the ob- 
served Q profiles, while the discrepancies for the somewhat 
uncertain slice 3 are less than 2a (Hovatta et al. 2012). The 
model has thus successfully reproduced the main features 
of the observed polarization. The total intensity fits were 
much less successful, with the model profiles generally being 
appreciably more symmetrical; differences between the ob- 
served and model / profiles were within 2a for the 6-cm Slice 
3, but frequently exceed 3a for the other slices considered. 
Whether this discrepancy between the model and observed 
I profiles is due to other factors giving rise to asymmetry in 
the observed I profiles or unsuitability of the model is not 
clear. For this reason, the results presented here should be 
considered tentative. 

Formally good model fits were not expected using the 
procedures described in this paper; the main goal was to 
determine whether a simple helical-field model such as that 
considered here is able to qualitative reproduce the shapes 
of the observed profiles, particularly in Q. In fact, there are 
many reasons why it would be unreasonable to expect a for- 
mally good fit to both I and Q. The structure of the jet may 
differ in many ways from the simple structure used here. For 



example, any deviation from perfect cylindrical symmetry, 
or an emissivity that changes with distance from the jet axis, 
would introduce features that could not be reproduced by 
the model. Some such physical deviations may be responsi- 
ble for the double-peaked I structure in the 3.6 cm profile in 
Fig. 6, which is poorly represented by the model. Attempt- 
ing to introduce further parameters into the model would 
result in fits that were much more poorly constrained by 
the data. Given the simplicity of the model used, its ability 
to reproduce the Q profiles so well is striking. 

Thus, the most important result of our profile fitting is 
that even the simple helical-field model considered here can 
reproduce most of the qualitative features of the polariza- 
tion profiles, lending support to the hypothesis that a helical 
magnetic-field component is present in the jet of Mrk 501 
on parsec scales. 

9.1 Curvature-Induced Polarization? 

Given the observed bends in the jet of Mrk 501, it is nat- 
ural to ask whether the observed polarization structure is 
associated with curvature of the jet. As the jet bends, the 
jet plasma is squeezed and stretched at the inside and out- 
side edges, respectively. If, as seems to be the case here, the 
magnetic field is highly disordered, but with a mean value 
of < Bf ong > greater than < B perp >, where Bi ong and 
B perp are the magnetic fields along and perpendicular to 
the jet, then squeezing the plasma along the jet axis will in- 
crease < B perp >, reducing the degree of polarization, while 
stretching the plasma in this same direction will have the op- 
posite effect. This effect would increase/decrease the degree 
of polarization at the outside/inside edges of jet bends. 

To investigate this, we constructed a map of the degree of 
polarization m at 6 cm, shown in Fig. 10. The observed frac- 
tional polarization is greatest at the edges of the jet (reach- 
ing values ~ 25 — 30%), as is expected for curvature- induced 
polarization. However, the locations of highest fractional po- 
larization are on the inside edges of local bends in the jet 
(the northern side of the jet between slices 1 and 2, and the 
southern side of the jet beyond slice 2), whereas curvature- 
induced polarization should be highest on the outside edges 
of these bends. Therefore, we find no evidence that the lon- 
gitudinal polarization observed at the edges of the Mrk 501 
jet is primarily induced by the curvature of the jet as it 
bends. 



10 CONCLUSIONS 

A method for comparing model and data profiles for total 
intensity and polarization in astrophysical jets has been de- 
scribed and demonstrated. The method has been used to 
compare observed profiles of the parsec-scale jet of Mrk 501 
with the predictions of a model in which a cylindrical jet is 
permeated by a magnetic field with a uniform helical com- 
ponent and a disordered component. 

The jet of Mrk 501 shows several characteristics that are 
consistent with a helical jet magnetic field, most strikingly 
the spine-sheath transverse polarization structure observed 
in the region of our Slice 2. The best model fits obtained cor- 
respond to pitch angles in the jet rest frame 7' ~ 40 — 50° 
and viewing angles in the jet rest frame 8' ~ 83°. These 
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fits describe the polarization structure well, though the to- 
tal intensity fits are generally poorer. However, given the ex- 
tremely simple nature of the helical-field model, for example 
its rigid cylindrical symmetry and uniformity of emission, its 
success in reproducing the general features of the Q profiles 
is noteworthy, and suggests that further comparisons, using 
higher resolution data and data for other sources, would be 
worthwhile. 

The best fit values of the line-of-sight angle are very sim- 
ilar for all the analyzed 4 cm, 6 cm, 13 cm and 18 cm slices 
across the Mrk 501 jet, within the estimated la uncertain- 
ties (apart from the fitted viewing angle for the somewhat 
more uncertain slice 3, which differs from the other values 
by 2 — 3a). This suggests that, as expected, the large appar- 
ent changes in jet direction are in fact very small bends that 
are greatly amplified by projection. The fitted pitch angle in- 
creases from the first to the second 6 cm slice, bringing about 
the observed transition in polarization structure (from con- 
figuration 1 to configuration 3 as described in Fig. 2). The 
fraction of energy in the disordered magnetic-field compo- 
nent seems to decrease with distance along the jet. 

Together with the tentative superluminal speed reported 
by Piner et al. (2009), the estimate for the viewing angle 
in the jet rest frame obtained through the profile fitting, 
8' ~ 83°, enables determination of the viewing angle and 
jet velocity in the obsever's frame, 5 « 15° and /3 « 0.96. 
Although these values are somewhat uncertain in the case 
of Mrk 501, since this superluminal velocity was determined 
for a weak feature whose position is only poorly defined, the 
joint analysis of transverse polarization profiles and appar- 
ent superluminal speeds provides a new tool for disentan- 
gling 8 and /3. The jet of Mrk 501 must be fairly close to 
the line of sight, yet carrying out this analysis for the typ- 
ically subluminal speeds observed in the Mrk 501 VLBI jet 
yields viewing angles of about 60°; this essentially demon- 
strates that these subluminal motions must represent pat- 
tern speeds in a much more highly relativistic flow, rather 
than highly relativistic motion viewed at a very small angle. 

Calculation of I and Q profiles for the helical-field model 
is greatly simplified if the spectral index a is assumed to 
be unity. Whilst this is almost certainly incorrect (observed 
values are usually ~ 0.5) the assumption that a = 1 has 
been shown to have very little impact on the profiles, and 
hence on the values of the best-fit model parameters. 

The results demonstrate that this method provides a new 
approach to studying the magnetic fields in parsec-scale jets. 
We are in the process of identifying other AGN jets that 
are well resolved and straight, display clear transverse po- 
larization structure, and possess reliable component speed 
measurements, which we hope will prove fruitful subjects of 
analyses similar to those carried out here. 
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